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Abstract 

The transverse momentum distribution of bosons at hadron colliders 
is well described by a parton-shower model for small p± values, but not for 
large ones. This article is an attempt to give a better description of the 
distribution by using corrections derived from the matrix-element formal- 
ism. The parton shower for qq' — > has been modified to resemble the 
matrix elements of qq — > gW and qg — > q'W at large p± values. Compar- 
isons between different approaches are presented at y/s = 1800 GeV. The 
results are also compared with experimental data from the DO collaboration 
at Fermilab. 



1 Introduction 



The Standard Model is the theory which lies at the heart of modern particle physics. 
In this theory we distinguish between two kinds of particles: matter particles and gauge 
bosons. (A boson is a particle of integer spin.) The Standard Model describes the matter 
particles (quarks and leptons) and their interactions. There are four fundamental forces 
in nature; interactions between particles take place via the exchange of the gauge bosons 
corresponding to these forces (see table below). 

Interaction Gauge bosons 

gravitation graviton (spin = 2) 

electromagnetic 7 (spin = 1) (1) 

weak W + , W~ Z° (spin = 1) 

strong gj, i — 1 . . . 8 (spin = 1). 

The gravitational force acts on all forms of energy, but is so weak that it can be 
disregarded in particle physics. The electromagnetic force acts on all electrical charged 
particles and is responsible for holding the electrons and the nuclei together in the atoms. 
It can create and annihilate photons. The weak interaction accounts, amongst other 
things, for the /5-decay of nuclei. In the Standard Model, the electromagnetic and the 
weak forces have been unified into one force: the electroweak interaction. The strong 
force is responsible for holding together the quarks inside the nucleons, as well as the 
protons and neutrons inside the nuclei. The theory for the colour (strong) interaction is 
called quantum chromodynamics (QCD)|l|. 

The quarks and leptons are the fundamental particles that matter is made of. They are 
all spin-1/2 fermions. To each matter particle there corresponds an antiparticle, having 
the same mass and spin as its partner, but with all other quantum numbers (e.g. electrical 
charge) having opposite values. 

The quarks are most significantly affected by the (unified) electroweak force and the 
strong force. There are six known 'flavours' of quarks, which are placed in doublets called 
'families' or 'generations'. Thus, the three generations of quarks are: 




The quarks in the upper row all have electric charge +(2/3)e, while those in the lower 

row have charge — (l/3)e. From the point of view of the strong interaction, each quark 

carries one of three possible QCD colour charges: r (red), g (green) or b (blue). 

All coloured particles are bound inside colourless 'hadrons'. From experiments we know 

that the hadrons are made up not only by quarks, but also by gluons. The concept of 

partons is thus introduced as a common name for the constituents of the hadrons (cf. the 

concept of the nucleon in nuclear physics). The hadrons are subdivided into 'baryons' 

(hadrons with half- integer spin) and 'mesons' (hadrons with integer spin). 

The leptons are unaffected by the strong interaction. They too are placed in doublets and 

the corresponding three generations are: 
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The electron e, the muon /x and the tau r all have electric charge — e, while the neutrinos 
v are electrically neutral. 

The electromagnetic and the gravitational forces have been known and studied for a 
long time. When it comes to the strong interaction, since the gluons are confined within 
hadrons, they can only be studied indirectly. The bosons of the weak interaction are 
special as they are so massive; they have been much studied since their discovery, in 
1983. Until recently, they could only be observed in hadron colliders, but now they have 
also been produced at LEP2. 

In this article we focus on the transverse momentum distribution of the W boson in pp 
collisions. This is related to QCD corrections to the basic electroweak process in which 
W's are produced. These studies are mainly of interest as a test of our understanding of 
the QCD, but they implicitly influence our confidence in the Standard Model. 
There are two complementary methods of describing the pj_ distribution of the W: the 
matrix-element approach and the parton-shower one. The matrix-element formalism 
has the advantage of being exact (to a given order). Unfortunately, the calculations 
become more and more complicated as one allows for more and more partons. In the 
parton-shower formalism, on the other hand, one can describe events with an arbitrary 
number of partons. This method is approximate, though, and it can only be trusted in 
the collinear limit (small p± region). 

In this Diploma Work we start from the parton-shower formalism and show that it can 
be extended to hold for large p±_ values as well, by adding corrections derived from the 
matrix-element formalism. 

This article is organized as follows. Section 2 includes a more detailed description 
of the matrix-element and the parton-shower formalisms, a presentation of the cross 
sections for W production and a description of the 'backwards evolution' method for 
reconstructing a parton shower. In Section 3 we present a comparison between the 
parton-shower and matrix-element differential cross sections, as well as a detailed 
description of the modeling of the parton shower according to the matrix elements. The 
results are presented in Section 4. Finally, Section 5 contains a summary of this article. 



2 Theory 

In the simulations of this Diploma Work, protons and antiprotons are allowed to collide 
at 1800 GeV in the CM reference frame, like it was done in the experiments at Fermilab. 
One can better understand why such high energies are advantageous to produce a W 
particle of mass ~ 80 GeV if one remembers that it is actually a quark and an antiquark 
- and not the whole proton and antiproton — that produce the W. The quark and the 
antiquark carry only a fraction of the energies of the proton and the antiproton. 

In QED it is well known that an electrically charged particle which is accelerating 
emits photons, for instance as is the case when an electron moves within a magnetic 
field. In general, high-energy interactions between charged particles are accompanied by 
photon emission. This phenomenon is known as bremsstrahlung. 

In a similar way, in QCD, scatterings between particles can give rise to new partons 
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being produced. Unfortunately ambiguities arise when attempts are made to distinguish 
between bremsstrahlung and hard-scattering partons. In a general process ab —>■ cd, 
partons can be radiated both by the incoming a and b (initial-state radiation), as well 
as by the outgoing c and d (final-state radiation) partons. Furthermore, interference 
between the initial- and the final-state radiation can arise, and these effects are not 
negligible. To simplify matters one can choose to consider processes where, for instance, 
only initial-state radiation is possible, such as qq' — > W + — > e + u c or qq — > Z° — ► e + e~. It 
is this initial-state radiation that is responsible for the transverse momentum spectra of 
the W and Z. (Actually, the shower initiators also possess some transverse momentum 
because they were confined within the protons. The effects that these so-called primordial 
p± have on the spectra are small, and will not be considered here.) 

In this article we study the transverse momentum distribution of the W bosons. The 
existence of the W particle was first proved in 1983 at CERN, where very energetic protons 
and antiprotons were allowed to collide. The W particle is unstable since it can decay into 
less massive states. Thus it has a characteristic lifetime r and, because of the Heisenberg 
uncertainty principle, an uncertainty in the mass. As a result, it is possible for W's of 
masses somewhat different from the nominal one to be produced. 



2.1 Mandelstam variables 

In a general process of the form A + B — > C + D the following invariants — so-called 
Mandelstam variables — can be defined 

s = (pa + Pb) 2 , (4) 
t = (PA-Pc) 2 , (5) 
w = (Pa-Pd) 2 , (6) 

where Pi is the four-momentum of particle i. The variable s is simply the square of the 
total energy in the CM system. If m» represents the mass of particle i then we obtain the 
relation 

s + t + u = Y / m 2 i . (7) 

i 

If one actually knows the subprocess a + b — > c + d, one can define the corresponding s, i 
and u by using p a , pb, p c and pd as arguments (see figure [I]). For instance if A is a proton 
and B an antiproton then a can be a quark and b an antiquark. 
One can define some dimensionless variables X( according to 

Pa = xtp A , (8) 
Pb = x 2 p B , (9) 

relating the incoming four-momenta of the process and the subprocess. This can be done if 
one assumes that a (b) moves along the direction of A (B), which is a good approximation 
at high energies. In order to derive a useful relation between s and s, consider the process 
pp — > WX with the corresponding subprocess qq' — > W. We have 

s = (Pp + Pp) 2 , (10) 

S = (Pq+Pq) 2 , (11) 
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Figure 1: The process A + B — > C + D. The circled region marks the corresponding 
subprocess a + b — > c + d. 



or 



s = 2mj + 2p p pp, ( 
s = (xip p + x 2 pp) 2 = {x\ + x 2 2 )m\ + 2x 1 x 2 p p Pp- ( 

By noticing that s > s and s ~ ^> m^, one obtains from equations flT2|) and ( |13"D 



(12) 
(13) 



s 



(14) 



Alternatively, it is possible to take the equation 



s = 



X1X2S 



(15) 



as a starting point and use it in the definition of x\ and x 2 - 
2.2 Matrix elements versus parton showers 

In order to describe interacting fields, in which particles can be scattered, created 
and annihilated, one has to solve the very difficult equations which arise. Considering 
in particular QED and QCD (which both have gauge bosons that are massless), one 
approach is to work within the framework of perturbation theory. Technically this means 
that the Hamiltonian of the system is divided into a free-field term plus an interaction 
term. For a sufficiently weak interaction, this last term can be treated as a perturbation. 

A matrix element expresses the probability amplitude for a scattering process 



Here \i) represents the initial state long before the scattering occurs, specifying a definite 
number of particles and their properties when they are far apart from each other. Similarly 
\f) describes the final state, after the particles have interacted and when they are far apart 
again. S is the so-called scattering matrix (S-matrix) JJ. 

Once the initial state and the final state have been fixed, the cross section is obtained by 
summing and then squaring all possible Feynman diagrams that contribute to the process, 
in the usual quantum mechanical manner. The more vertices we include in the diagrams, 



M = (f\S\i). 



(16) 



the more difficult the calculations become. 

To calculate the contribution from each diagram one makes use of the Feynman rules. A 
factor of ^/a em (^/aQ corresponds to each QED (QCD) three-particle vertex. The four- 
gluon vertex of QCD corresponds to a a s factor and is thus more suppressed. 
The cross section can be expanded in a a em (a s ) series. In the case of QED, the expansion 
terms decrease at a relatively fast rate. Consequently, it is a good approximation to keep 
only the first few terms of the expansion, as the rest terms are negligible. For QCD, the 
terms do not decrease as fast, and more terms in the expansion are of significant size. 
In QED the cross section is expanded in terms of the small dimensionless fine-structure 
constant a em , which measures the strength of the coupling between the electron and the 
photon 

e 2 1 

a em = — ~ . (17) 

4ne hc 137 V ; 

For QCD one can take a similar approach by expanding in the corresponding so called 
strong-coupling constant. This in not really constant, but is 'running' with the energy. 
(At this point we should add that in QED, a em is not a constant either; the expression 
given above is the value at Q 2 = 0.) To first order, it is given by 

a s {Q . 



(33-2n f )ln(Q 2 /AQC D )' 

where Q is the energy scale for the process, rif is the number of quark flavours available 
at the actual energy (usually 4-5) and Aqcd ~ 200 MeV is the QCD scale parameter ||. 
The matrix-element description has the advantage of being exact (to a given order). It 
can be used successfully for processes where the outgoing particles are well separated. 
Unfortunately it is not so easy to apply in the collinear limit, since a s becomes larger 
and the series converges slower. Another drawback of the matrix-element description 
is that the calculations become complicated very fast as one allows for more photons 
or gluons in the processes. Also, there are uncertainties in the choice of the Q 2 scale for a s . 

A complementary approach to the matrix-element description is the parton-shower 
one. The building blocks of the parton shower are branchings of the form a — > be. These 
can be repeated, forming a tree-like structure. Each branching vertex is associated with 
some relative transverse and longitudinal momentum between the partons b and c. For 
an a — > be branching, we let the variable z represent the fraction of the longitudinal 
momentum of the incoming particle a that particle b is taking. 

The three basic QCD branchings q — > qg, g — > gg and g — > qq are shown in figure ^| (the 
branchings in the QED case are analogous, except that there is no three-photon vertex, 
since photons do not interact with each other). The branching probability P a _+i, c (z) for 






q -> qg g -> gg g -> qq 



Figure 2: The basic QCD branchings 
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each case is given by the so-called Altarelli-Parisi (AP) splitting kernels 



4 1 4- z 2 

*W*) = 3^7. (19) 

W*) = 3 (1 ; ( f_"^ ))2 , (20) 

P g ^(z) = n^(z 2 + (l-z) 2 ). (21) 

These can be derived from the expressions of 3-parton-matrix-elements cross sections, 



by taking the limit of two partons being collinear. We will return to this in Section |3T 
Notice that the second equation corresponds to two identical particles in the final state. 
The factor 3 therefore has to be replaced by 6 in the case that one is asking for the 
probability of finding a gluon in the final state, given that its mother was a gluon. The 
first expression is singular in the limit z — *■ 1, which corresponds to a 'soft' (of low 
energy) gluon being emitted. Similarly, the second expression is singular both as z — > 1 
and as z — > corresponding to any one of the two emitted gluons being 'soft'. 
In the described model of the parton shower, the different branchings are regarded as 
independent. As a consequence the total amount of evolution is overestimated. In the 
case of initial-state radiation, coherence effects can be taken into account by demanding 
that the virtualities associated with the partons on the 'main chain' - the chain that 
starts with the initiating parton and ends with the scattered one - are ordered, with the 
largest being closest to the hard scattering (cf. Section |2.4|) . 

The parton-shower description has the advantage that one can allow for an arbitrary 
number of particles both in the initial and the final state. Also, it gives a good description 
in the collinear limit. Unfortunately the description is approximate. 

There are two alternative approaches one can adopt in describing the change in the 
number of partons, as a function of the resolution. The first picture is an exclusive one: we 
follow the original parton as it repeatedly branches. We will return to this in Section E^. 



The second approach is the parton-density picture: we give an inclusive description of the 
number of partons of a certain kind. Consider the proton, for instance. Nonperturbatively, 
it is a bound state and its parton composition is not known beforehand. Given some 
(fitted) distribution at a low Q 2 scale, the change in the number of partons b is given by 
the so-called DGLAP evolution equations 

^-^H^^"^- (22) 

This inclusive picture of showers is also used in the matrix-element formalism, but, unlike 
in the normal showers, no is assigned to the partons. If we put t = ln(Q 2 /Aq CD ), the 
above equation can be rewritten as 

^ =£ / 1 ^^) Sa (J, t)P ^ Uz) . (23 ) 

QI „ Jx X ATI 



The functions fi(x, t) are the so-called parton distributions, which measure the probability 
of finding a parton i carrying a fraction x of the total momentum of the incoming hadron. 
The variable Q 2 represents the scale for the process: the higher the Q 2 value, the finer 
a structure in the hadron one can distinguish. Here x' and zx' = x are the momentum 
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Figure 3: The sharing of the momentum fractions x', x and z. 

fractions that the partons a and b, respectively, take from the total hadron momentum, 
with z defined as before (see figure ||). 

To better understand the meaning of the DGLAP evolution equation, consider the case 



when b represents a quark q. Then equation (|23D becomes 



The first term in ( p4|) expresses the fact that the quark with momentum fraction x could 
have originated from a parent quark with a larger momentum fraction x' which has 
emitted a gluon. The probability for the emission is proportional to a s P q ^ qg (,2 = 4) and 
one has to integrate over all possible momentum fractions x' > x of the parent. The 
second term can be interpreted in a similar way, but now the parent is a gluon. 

To summarize, we conclude that the matrix-element description is useful for the de- 
scription of the p± spectrum at large values (when a s is small), while the parton-shower 
description is a good candidate in the small p± region (collinear limit). 

2.3 Cross sections for W production 

In a collider, the number of events of a particular kind is given by 

N = aL (25) 

where the integrated luminousity L is a characteristic of the accelerator and the total 
cross section a is a property of the process in question. 

If protons and antiprotons are allowed to collide at the CM energy y/s, then the lowest- 
order process that produces W's is qq' — > W (see figure f|). The total cross section is 
given by the expression 

c = o"o / / dxi dx2 F(xi,X2) 5 ( X1X2 — ) , (26) 




Figure 4: qq' — > W. 
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where 



7r 2 a P 



a ° ~ 3 sin 2 fl w mfv' (27) 
^(xi,x 2 ) = £{*i/*(*i,QV./?(z2,Q 2 ) +x 1 /P(x 1 ,g 2 )x 2 /P(x 2 ,g 2 )}|Vj| M | 2 .(28) 

qq' 

V is the Cabibbo-Kobayashi-Maskawa (CKM) matrix. It is a 3 x 3 unitary matrix 
containing three independent mixing angles and one phase. The angles are generalizations 
of the Cabibbo angle, describing the possible mixings between the (u,d), (c,s) and (t,b) 
doublets. As an example, one would naively expect a process like ud — > W to be allowed 
but us — > W to be forbidden. Because of the mixing, however, even the latter process is 
possible, although at a reduced rate. 

The mass-distribution spectrum for the W has a Breit-Wigner shape, with a peak at the 
nominal mass rs 80 GeV. The Breit-Wigner distribution is given by 

mx7Y-\ 1 m w r w 

Bff( ''°; (i-o-Ki', ' (29) 

where T = 1/r is the full width at half maximum. It expresses the fact that there is a 
non-zero probability of producing W's with the 'wrong' mass. 

Unless one is explicitly interested in the W mass distribution, it is a good approximation 



to replace the actual Breit-Wigner resonance by a 5-function, as we did in equation fl26|) 



In the collinear limit, processes like qq' — > gW and qg — > q'W are already included 
in the expression for the cross section given in equation (p6|). (The Q 2 dependence of 
the parton distributions fi(x,t) from equation implicitly means that the x scaling 
violations are included.) Here we look at the corresponding matrix elements, relevant at 
large p± in particular. For the qq' — > gW process, the contributing Feynman diagrams 
are shown in figure |5|. These diagrams interfere and the cross section is given by 




Figure 5: qq' — > gW 



a = [[[dx^di ^^ (30) 

J J J X\%2 dt 

da 4« s mfv i 2 + u 2 + Ini^s . 

Wh6re Yt = ao 3 271^ Tu ' (31) 



and where a and F(xi,x 2 ) are as in equations (p7|) and (p8|). Notice that da/dt is 
symmetric in t and u. 
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s-channel graph 



w-channel graph 



Figure 6: qg — > q W 

In the case of qg — > q W the contributing diagrams are shown in figure []. The cross 
section is formally given by the same expression as for the qq' — > gW case (equation 
(PD|)). Because we here also have a gluon in the initial state, the terms in the expression 
for F(xi,£ 2 ) will include structure functions of the type f^(x,Q 2 ). Also, da/dt is here 
given by 

da_ 1 Os raw s 2 + it 2 + 2m w t 

dt ~ a ° 2 27r P ■ W 

Notice that da/dt is now symmetric in — s and u. (With a different notation for the 
particles in figure ^ i.e. if particles 3 and 4 are interchanged, the t and u would have 
been interchanged in the expression da/dt above). 

The W + can decay through one of the following main channels: 

W + — > e + v e , ud, cs, r + z/ r , tb, (33) 

and the W~ decays in a similar way. (Actually, more channels are possible because of 
the mixing effect mentioned above.) In reality the tb channel is not allowed, since W is 
not heavy enough. Because quarks can have three different colours, one has to remember 
that a decay like W + — > ud actually corresponds to three different channels W + — ► u r d r , 
W + — > Ugd g or W + — > Ufedb. 

The probability of W to decay to a final state i is given by T l w /T^, where 

rW = ^^^(N c ) l |Vf KM | 2 m w . (34) 



12 sin^ 6 



w 



(N c )j is a colour factor which equals 3 for quarks and 1 for leptons. r w * is the total width 
for the W decay 

r w t = ^r w ^2GeV. (35) 

i 

Thus, one obtains for the lifetime of the W 

T = h/T ~ l(T 25 s. (36) 

2.4 Backwards Evolution 

In the QCD shower description, the process qq' — > W can look schematically as in 
figure [f[ Parton showers are allowed in both the initial and final state. The hard 
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process (highest virtuality) is the part of the diagram where the W is produced; it is the 
propagator that sets the scale for the process. The maximum virtuality is set equal to 
the square of the mass of the propagator, Q max = m 2 ^. Both in the initial- and final-state 
cascade the virtuality increases towards the hard process. 

In this article we are interested in the transverse momentum for the W's produced, so 
only the initial-state parton shower is relevant. 
In describing the initial-state radiation, one speaks of a 'main chain' of consecutive 

initial-state parton shower final-state parton shower 




Qo *" Qma,x — m w " m o 



Figure 7: Initial-state and final-state parton showers. In both cases the virtualities are 
increasing towards the hard scattering. 

branchings, which connects the initiating parton to the scattered one. In the successive 
branchings a — > be of the cascade, the partons a and b are on the main chain. Conven- 
tionally, these partons are assigned a spacelike virtuality, while c is assumed to be on 
mass shell (or of timelike virtuality, in which case it will initiate a final-state cascade). 
Since the (spacelike) virtuality of the cascade is carried by one single parton at a time, 
it is possible to equate this with the Q 2 of the cascade, used e.g. in the DGLAP equations. 

If the cascade is evolved in the 'forward' direction, starting with the cascade initiators 
and ending with the hard scattering, it is not beforehand possible to decide which partons 
are spacelike, so this and other complications arise. The problem is solved by adopting 
the 'backward evolution' scheme Q], where we start by choosing the hardest process 
and work our way back to the shower- initiating partons. The cascade is reconstructed, 
by making use of the AP splitting kernels and the Q-dependent structure functions 
(cf. Section |2.2| ). During the evolution, one proceeds to lower and lower values for 
the virtuality. A minimum limit for the Q 2 is set so that the evolution will stop when 
Q 2 = Ql where Q 2 ~ (lGeV) 2 §. 

The starting point for the backwards evolution is the DGLAP equation (cf. equa- 
tion (H)) 

^=E/Tf*1W.), (37) 
where, as before, the variable t = ln(Q 2 / Aq CD ) . Now we are working backwards in time, 
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away from the hard scattering. Therefore the DGLAP equation expresses the rate at 
which partons b of momentum x = zx' are 'unresolved' into partons a of momentum x', 
as we take a step dt backwards. 

The corresponding relative probability is dP^/dt = (1//&) (dfo/dt). The cumulative effects 
of many small dt steps are summarized in the so-called Sudakov form factor, which gives 
the probability that parton b stays at x from t max to t < t max : 



The differential probability distribution dSb/dt that a branching a —>■ be actually takes 
place between t and t — dt is 



Once the Sudakov form factor is known, the parton shower can be reconstructed. At 
each branching, the values for t (also defining the virtuality of parton b), z (the splitting 
variable) and a (the flavour of the mother parton) need to be found. The variable t is 
chosen according to the differential probability distribution dSb/dt. For a given t, a is 
chosen according to the branching probabilities — the x' integrals. Finally, with t and a 
known, the probability distribution for z is given by the mentioned integrand. 
The process of reconstructing the cascade is iterative; when going from one branching 
a — > be to the next a' — > b'c' (where a = b'), the actual t value of the parton b is taken as 
the upper limit t max for the parton b' . Thus, the virtualities of the partons on the main 
chain decrease during the evolution. The opposite is true for the momentum fractions x. 
(Remember that we are evolving backwards.) 

Finally, some words about the interpretation of the variables z and s. We first consider 
the case of collinear and massless particles and use the notation of figure |8|. For the partons 



Figure 8: Schematic picture of initial-state parton showers. Partons 1 and 2 take part 
in the hard scattering. 

1 and 2 that form the W we have as before 




(38) 




(39) 




(Pi +P2) 2 



m w(= Sq). 



(40) 
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On the other hand, if parton 1 is unresolved, we have that 



(P3+P2) 2 = S. (41) 

The splitting variable z relates the four-momenta of parton 1 and its mother, parton 3, 
according to 

Pi = zp 3 . (42) 
By multiplying both sides with p%, this can be rewritten as 

(P1+P2) 2 = HP3+P2) 2 - (43) 

(Here we have used the assumption that the particles are massless.) From this and 
equations (HOI) and ( |4T| ) we obtain 

z = (44) 
s 

This expression is used iteratively in the 'backwards-evolution', thus also for non-collinear 
emission. We will use this expression as a definition for z, when we will make translations 
between the parton-shower variables z and Q 2 and the matrix-elements ones s, t and u 
(cf. Section [3|). 



3 Modeling the parton shower according to the ma- 
trix element 

The parton-shower description has the advantage that one can allow for an arbitrary 
number of particles. This is in contrast with the matrix-element description, where the 
calculations become complicated very fast as the number of particles increases. Also, we 
expect the parton-shower formalism to give a good description of the distribution at low 
p± values, and the matrix-element formalism to be useful in the high p± region. 
The main goal of this work is to take the lowest-order matrix element of qq' — > W 
with initial-state parton shower, and modify it so that it can be used as an alternative 
description to the higher-order matrix-element at large p± values as well. 

From a QCD point of view, the lowest-order matrix element qq' — > W is of zeroth 
order in a 8 . Since we want to describe the W transverse momentum, this 'naked' process 
is uninteresting. Instead we have to allow initial-state radiation. In this Diploma Work, 
we start by considering the qq' — ► W with initial-state parton shower, and introduce 
corrections derived from the matrix elements of the processes qq' — > gW and qg — ■> q'W. 
We show that the thus modified parton-shower description can be extended to hold for 
large p± values as well. In the rest of this article, the following notation will be used: 

• PS : qq' — > W with initial-state parton shower, 

• ME : qq' — > gW and qg — > q'W. 

Experimentally, one might wish to consider the two processes qq' — > gW and qg — > q'W 
together, as they are both of first order in a s . We choose to make distinctions between 
them, and we will see in Section [| that this separate treatment is well motivated indeed. 
Thus, for comparison reasons, we divide the PS into two parts: 
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Figure 9: qq' — > W with initial-state parton shower. The circled region is similar to the 
matrix element of qq' — > gW. 



Figure 10: qq' — > W with initial-state parton shower. The circled region is similar to the 
matrix element of qg — > q'W. 

1. PS similar to qq' — > gW. 

The two Feynman diagrams that contribute to this ME are shown in figure [|. For 
the PS to be similar to the above ME, the most virtual quark (closest to the hard 
scattering) must have come from a q — > qg vertex, as is shown in figure []. 

2. PS similar to qg — > q'W. 

The two Feynman diagrams that contribute to this ME are shown in figure For 
the PS to be similar to the above ME, the most virtual quark (closest to the hard 
scattering) must have come from a g — > qq vertex, as is shown in figure [K]. 

3.1 Comparing the ME and PS differential cross sections 

In this section we compare the differential cross sections for the PS and the ME. An 
expression for the PS differential cross section is derived from the ME one, as illustrated 
below. 

1. qq' — > gW compared to the corresponding part of the PS. 

The Feynman diagrams for this ME are shown in figure |^. The total cross section 
is given by equation (|30|) where 
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In order to derive the corresponding differential cross section for the PS, we first 
rewrite Ame, by replacing the matrix-element variables s, t and u with the parton- 
shower variables z and Q. 

The t-channel graph in figure |5| is similar to the parton-shower diagram (if only the 
hardest gluon of the cascade is considered) shown to the left in figure §. We can 
make the following 'translations' 



s = (Pi+P 2 ) 2 = ^, (46) 

i = (Pi-P3) 2 = -Q 2 , (47) 

I — z 

u = m\j - s -t = Q 2 rn^ff. (48) 



In the first equation we have used the expression for z of equation (|44]) . In the 
second one we have equated t with the virtuality of the cascade. In the third we 
have used equation (|7]). 

The expression for A M e can now be rewritten as 



A 



4 t 2 + u 2 + 2m 



ME 



3 tii 

4 - 2 l -^Q 2 m 2 w + 2Q 4 
3 g 2( i^ m 2 v _g 2) 

g 2 3 i - z z ■ 1 ] 

In the last row we have assumed small virtualities Q 2 -C m^, which correspond to 
the region where PS can be trusted. As expected, we recover the AP splitting kernel 
P q ^ q g(z) and also the 1/g 2 factor present in the DGLAP evolution equation (cf. 
equation (^)). Thus we see that the parton-shower-formalism expression is indeed 
recovered at small p±. (The extra factors present in the expression cancel when the 
kinematics are considered and the integration of equation (|3~0|) is performed.) 
This PS expression for A, which we will call Apsi, can again be expressed in terms 
of s, t and u by making use of equations ([46]), ( [4TD and (|48|). We thus obtain 



o (t + u)t 

The w-channel graph in figure [5] is similar to the parton shower diagram (only the 
hardest gluon of the cascade is considered) shown to the right in figure ^|. The 
corresponding 'translation' equations are obtained as above, simply interchanging t 
and u. Thus the expression Aps 2 is 

_ 4 s 2 + r< 

Ap S 2 - o a , rv^ • y bi > 
o [t + u)u 

Since the parton-shower description is additive (no interference), the total parton- 
shower answer is simply A PS = A PS1 + A PS2 , i.e. 
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(52) 



Finally, we obtain an expression for the ratio R(s, t) between the ME and the PS 
differential cross sections: 



R(s, t) 



(da/dt) 
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One can show that 
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< fl(s,t) < 1. 



(53) 



(54) 



2. qg — > q'W compared to the corresponding part of the PS. 

The calculations are similar to the previous case. The Feynman diagrams for this 
ME are shown in figure [| The total cross section is given by equation (|30|) where 
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(55) 



As before, Ame is to be rewritten by replacing the matrix-element variables s, t and 
u with the parton-shower variables z and Q. 

The w-channel graph in figure |6| is similar to the parton-shower diagram (if only 
the hardest gluon of the cascade is considered) shown in figure [10]. (The s-channel 
graph in figure ^| does not have a parton-shower correspondent.) The 'translation' 

T7D and 



equations are similar to equations (f46|) , ( |47D and (f48|) , but with the roles of t and u 
interchanged. Expressed in the new variables, Ame becomes 
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(56) 



In the last row we have again assumed small virtualities Q 2 -C m^, which corre- 
spond to the region where PS can be trusted. As expected in this case, we recover 



the AP splitting kernel P, 



qq v 



and also the 1/Q 2 factor present in the DGLAP 



evolution equation (cf. equation (|22|) ). This PS expression for A is rewritten in terms 
of s, t and u 
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The parton-shower answer in this case is 
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The expression for the ratio R(s, t) between the ME and the PS differential cross 
sections is thus given by 

R , g * = (dff/d£) M E = s 2 + u 2 + 2m^i ^ 
(d(j/dt) PS s 2 + 2m^« + 2my" 

Again, one can show that 

1£fl(i - i)£ 5ivfe) <3 ' (60) 

One possible motivation to the fact that R(s, t) is bigger in this case is that the 
s-channel graph is absent from the PS description. 

3.2 Description of the modeling 

The W transverse momentum distribution is described well by the lowest-order matrix 
element with parton shower in the region of small transverse momenta, but the descrip- 
tion is poor for large ones. This is closely related to the fact that, conventionally, the 
maximum value for the virtuality associated with the 'main chain' in the shower, should 
be of the order of the hard-scattering scale. In the case where a heavy resonance particle 
is formed, it is natural to let its mass set the scale for the process. Thus the maximum 
virtuality is set equal to the square of the mass of the resonance particle, in our case: 
Qmax = m w- This forces the distribution to fall abruptly for p± ~ Q max = m w . 

The main goal of this Diploma Work is to modify the PS description of the transverse 
momentum distribution, so that it can be used as an alternative approach to the higher- 
order matrix-elements description, also for large p±. The modeling is done by adding 
corrections derived from the ME. 

The events are simulated by using the event generator PYTHIA, written in FORTRAN77. 
In the generator, the initial-state parton shower is implemented as a Monte Carlo simula- 



tion. This makes use of the 'backwards evolution' scheme (cf. Section |2.4|) . The changes 
that we made have been implemented as a subroutine, which is a modification of the 
already existing subroutine PYSSPA. 
The PS is modified in two steps: 

1. The maximum virtuality is increased. 

As mentioned above, in the case of a resonance particle being formed, one conven- 
tionally sets Qmax = m w- This is also the default value that is already implemented 
in PYTHIA. 

On the other hand, when two hadrons are allowed to collide at a cm. energy of 
Eqm, the kinematic limit for the virtuality of the most virtual parton on the 'main 
chain', which actually takes part in the hard scattering, is of the order Eq M . We 
choose to increase the upper limit for the virtuality, by setting Q^ ax = (-Ecm/2) 2 . 

2. ME corrections are introduced. 

The differential probability distribution dSb/dt of equation ( pEOP is modified, by using 
the ratio R(s, t) between the ME and the PS differential cross sections, derived in 
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Section |Q| (recall that the ratio R(s, t) can be rewritten as a function of the parton- 
shower variables z and t). We modify as follows 



dt a J x X ' 27T fb(x,t) 

x exp (- I dt £ I — ^ hM P^ bc {z) R(s, t) j . (61) 

(Remember that the expressions for R(s, t) differ for the two different parts of the 
PS.) This is done only for the branchings closest to the hard scattering, on either 
side (cf. figures [5] and [TD]) . The evolution variable t is chosen according to the 
corresponding probability distribution §b(x, t ma3l , t), by solving the equation 

/•tmax dSb 

Sb(x,t max ,t) = j — = R, (62) 

as usual in a Monte Carlo approach. Here R is a random number between zero and 
one, chosen from a uniform distribution. To solve equation ( |62"D we make use of the 
veto algorithm 0: the (exact) expression for R(s, t) is replaced by R max (s, t) (where 
R(s,t) < Rmax(s,t)) and t is selected; thereafter the ratio (R(s,i))/(R ma , x (s : i)) is 
used to decide whether to keep this value for t, or to evolve further. In the case 
of the PS part that is similar to qq' — > gW, R(s, t) is given by equation ([53]) and 



-Rmax('M) = 1 (cf- equation (0)). Analogously, in the case of the PS part that 
is similar to qg — ► q'W, R(s,t) is given by equation ([59|) and -R max (s, £) = 3 (cf. 
equation (§H)). 

In addition to the final runs, intended to show the full complexity of the PS and 
allow a direct comparison with data, a number of test runs are performed to check the 
implementation of the ME corrections to the PS. 

In order to check the kinematics of the as-above modified PS, in the framework of the 
correction factors, we make the PS as similar as possible to the ME. Thus we artificially 
change the modified PS: we strongly suppress the probabilities of emitting all partons of 
the shower, except the most virtual (closest to the hard scattering) parton (cf. figure ^| 
and figure |10|). (For technical reasons, we can not generate events that are guaranteed to 



have only one parton in the shower; in our treatment we thus suppress the rate of multiple 
emissions to very close to zero, and neglect all additional partons.) 

To isolate the important features, that we want to study and compare for the PS and ME, 
we choose to generate events in which (1) no associated timelike showers are allowed, i.e. 
the emitted partons in the initial-state parton shower are put on the mass shell and (2) 
no final-state radiation is allowed. 

Finally, we turn to the changes introduced for technical reasons. The main problem, when 
trying to compare the PS with the ME transverse momentum spectrum for the W, is that 
of too little statistics at large p±. To find a way out of this problem, we first generate 
events in the usual way (with the above-mentioned changes) to cover low p± values, and 
then concentrate on the large p± region. In the ME case, we generate as above, but 
impose a minimal p± value on the generated processes. In the PS case, the situation is 
somewhat more complicated, and we need to introduce an artificial trick: we generate as 
above, with the additional requirement that the probability of emitting the most virtual 
(closest to the hard scattering) parton of the shower is strongly enhanced. Finally, these 
spectra are artificially connected at an intermediate p± value for the PS and the ME, 
respectively. 
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4 Results 



The main topic of this work is the study of the transverse momentum distribution of 
the W bosons. To do this we generate pp collisions at ^fs = 1800 GeV. Experiments 
involving pp collisions at y/s = 1800 GeV have been conducted at the Fermilab Tevatron 
Collider. 



We begin by presenting the test runs, which are performed in order to check the imple- 
mentation of the ME corrections to the PS. (To do this, the artificial changes mentioned 
in Section |T2] are also introduced.) 

In figure |11] the (unseparated) PS is compared with the ME. As expected, the distribution 
of the unchanged PS (short-dashed curve) falls abruptly at p± ~ Q m &x = mw- Notice 
that the increase of the maximum virtuality to Q^ax = (-^cm/2) 2 (long-dashed curve) is 
enough to prevent this kind of behaviour. Even more, the agreement with the ME (the 
dotted curve) is indeed very good, already by this simple change. The difference between 
the ME and the two PS curves, in the low p± region, is of technical, rather than physical, 
nature. (The ME has been generated with a p±_ m i n = 10 GeV.) 

Figures |T2| and [□] show the PS similar to qq' — > gW and to qg — > q'W, respectively, 
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Figure 11: PS (unseparated) compared to the ME. 



each compared with the corresponding ME (the dotted curves). Just as in the case of the 
unseparated PS, the distributions of the original PS (short-dashed curves) drop abruptly 
at p± ~ Qmax = m w- After increasing the maximum virtuality (long-dashed curve), the 



PS of figure [12] lies above the ME, while the opposite is true for the PS in figure |13|. 
(These effects seem to cancel each other out and are therefore not noticeable in figure [IT].) 
Clearly, the two parts of the PS need to be treated separately. The corrections derived 
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Figure 12: PS similar to qq' — > gW compared to the corresponding ME. 
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Figure 13: PS similar to qg — > q'W is compared to the corresponding ME. 



19 



from the corresponding ME in each case are necessary to give a better matching of the 
PS to the ME. The full curves include also these corrections. These curves agree quite 
well with the ME for all p± values. In the large p± region, though, they seem to lie sys- 
tematically slightly below the ME. For reasons of limited time, we have not yet examined 
this closer. 

Finally, we consider the behaviour of the correction factors R(s,t). From figure |14| we 
conclude that R(s, t) is close to unity for most of the events (notice the logarithmic scale), 
but there are indeed tails of the distributions. Also, the R(s, t) values for the two parts 
of the PS lie in the expected intervals (cf. equations (f>|) and flSOf)). As can be seen 
in figure [H| the R(s,t) — > 1 as p± — > 0, in accordance to what one expects, since it 
corresponds to the collinear limit. This is also the region with most branchings, which 
explains the peaked behaviour of the R(s, t) distribution. There are variations of the 
R(s,t) for p± > 0, but eventually, for p± £ mw, the mean value of the R(s,t) becomes 
fairly constant; this is true for both parts of the PS. 



We now present the full PS, in all its complexity (the full PS includes the two-steps mod- 
ification described in Section |3.2| , without any additional changes, unless specified). 
In figure [16] the full PS (long-dashed curve) is presented together the experimental data 
from the DO collaboration (the error bars account for both statistical and systematic 
uncertainties). The agreement is relatively good, for large p±, but for pj_ below ~ 20 
GeV, the PS is slightly shifted to lower p± values, compared to the data. To quantify this 
effect, we also perform runs with a larger primordial p± than the default value of 0.44 
GeV. The distribution corresponding to primordial p± = 4 GeV agrees very well with the 
data. Notice that the increase of the primordial p± mainly affects the distribution in the 
low p± region, leaving it fairly unchanged in the high p± region. 

In figure [L7] we present the full PS (short-dashed curve) together with the ME (continuous 
curve). The distributions agree in the small p± region; at large p± we can not really draw 
any conclusion, since the full PS alternates in being above and below the ME. 
Finally, in figure 18, we compare the full PS with the original one. The plot represents 
the difference between the full PS and the original PS, normalized to their sum. With the 
exception of a small region close to zero (p± ^ 5 GeV), the full PS distribution lies below 
the original one, up to p± ~ 40 GeV. For larger p±, the full PS lies above the original PS. 



5 Summary and Outlook 

In summary, we conclude that the parton-shower formalism can indeed be extended 
to hold in the large p± region as well, by introducing corrections derived from the 
matrix-element formalism. The parton shower for qq' — > W has been modified to 
resemble the matrix elements of qq — > gW and qg — > q'W at large p± values. During 
the modeling, the separation of the PS into two parts (each corresponding to a matrix 
element) was necessary. All pp collisions have been performed at yfs = 1800 GeV. When 
comparing the modified PS with experimental data from the DO collaboration, we found 
that the agreement was very good. 

The parton-shower formalism, modified as described in this article, can be applied in the 
case of the Z° bosons as well. 

A possible continuation to this work is to study pp collisions at ^ = 14 TeV. Such 
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Figure 14: The R(s, i) distributions for the PS similar to qq' — > gW and qg — > q'W, 
respectively. 
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Figure 15: The {R(s, t)) as a function of the W p± for the PS similar to qq' 
and qg — > q'W, respectively. 



gW 



21 



0.1 



0.01 



0.001 



0.0001 



1e-05 



1e-06 



PS, new 

PS, new, prim. pT = 4 GeV 
'9803003data' 



20 



40 



60 



80 



100 



transverse momentum for W (GeV) 



Figure 16: The full PS together with data from the DO collaboration. 
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Figure 17: The full PS and the ME. The original PS is also displayed as reference. 
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Figure 18: The difference between the full PS and the original PS, normalized to 
their sum, as a function of the W p±. 

collisions will be possible to perform at the LHC at CERN. For the pp runs at y/s = 14 
TeV that we have performed (not presented in this article), the discrepancies between the 
(separated) PS and the corresponding ME are larger than in the case of the pp collisions 
at y/s = 1800 GeV. We have not yet come so far as to be able to draw any conclusion 
on whether the differences are due to a kinematic mismatch, or some other physically 
nontrivial reasons. Once this is overcome, the approach offers hope to model p± spectra 
for signal and background processes with good accuracy by rather simple means, so it is 
of interest to study further. 
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